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tilt angle, collector and storage fluid stream velocities, 

and collector to storage heat exchanger dimensions. The 
procedure includes an accounting for economic parameters as 

an intimate part of the design process. The resulting 
methodology has been used for the design of solar energy 
systems which would use shelf item collectors for the purposes 
of determining the optimum design variable vector for a given 
Situation. The methodology could also be used on a limited 
basis for collector design optimization by exploring the 
effects of changing selected collector parameters on system 
performance. The methodology is coded in the FORTRAN computer 
language under the name SOLOAD-1 (SOLAR ENERGY OPTIMAZATION 
ANALYSIS OR DESIGN). ; 


Initial system trials indicate complete stability with 
minimal constraint activations. Based on the results of 
approximately fifty design experiments using SOLOAD-1, new 
findings concerning optimum collector tilt angle and an 
incariant optimum collector flow factor are suggested. 
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ABSTRACT 


‘A methodology for solar energy system design improvement has been 
developed and coupled with a constrained function optimization code 
resulting in an automated solar energy system design procedure. The 
scope of the methodology is limited to systems using flat plate 
collectors and water as the working fluid. 

Eight parameters have been included as independent design variables. 
The design variables included collector area, collector tile angle, 
collector and storage fluid stream velocities, and collector to storage 
heat exchanger dimensions. The procedure includes an accounting for 
economic parameters as an intimate part of the design process. The 
resulting methodology has been used for the design of solar eneray 
systems which would use shelf item collectors for the purposes of 
determining the optimum design variable vector for a given situation. 
The methodology could also be used on a limited basis for collector design 
Optimization by exploring the effects of changing selected collector 
parameters on system performance. The methodology is coded in the 
FORTRAN computer language under the name SOLOAD-1 (SOLAR ENERGY 
OPTIMIZATION ANALYSIS OR DESIGN). <— 

Initial system trials indicate complete stability with minimal 
constraint activations. Based on the results of approximately fifty 
design experiments using SOLOAD-1, new findings concerning optimm 
collector tilt angle and an invariant optimm collector flow factor 


are suggested. 
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NOMENCLATURE 


English Letter Symbols 


load heat transfer surface area, £t* 
collector area, fe? 

capacity, Btu/hr F 

collector cost, s/fte 

fuel cost, $/Btu 


systgm initial cost per unit collector area, 


$/ftte 
system ins$allation cost per unit collector 
area, $/ft< 


system lifecycle cost, $ 


system annual pumping power cost per unit 
collector area, S/fts yr 


system Lifecycle cost per unit collector 
area, $/fts 
storage tank cost per unit collector area, $/fte 


lifecycle operation and mgintenance cost per 


unit collector area, $/fte 


collector to storage heat _exchanger cost per 
unit collector area, $/ft2 


specific heat, Bty/lb F 


monthly diffuse solar insolation on a horizontal 
surface, Btu/ft2 mon 


collector loop inner diameter, ft 
collector loop out diameter, ft 
heat exchanger annulus inner diameter, ft 


value of the solar energy, $ 


th a ail call 


Pout 
D 


storage 


rh 


lifecycle cost, $ 


net savings, $ 


annual fraction of the energy load provided 
by solar energy 


monthly fraction of the energy load provided 
by solar energy 


collector efficiency factor 

collector flow factor 

system maintenance cost factor 

collector heat removal factor 

collector to heat exchanger efficiency factor 
present worth factor 

objective function 


collector loge flow rate per unit collector 
area, gpm/ft< 
inequality constraint 


sOlar insolation on a horizontal surface, 
Btu/hr ft 


extraterrestrial solar insolation ona 
horizontal surface, Btu/hr ft 


solar insolation ona gurface tilted towards 
the equator, Btu/hr ft* 


control flag for COPES 
heat exchanger length 

mass flow rate, lb/hr 
number of constraints 
number of design variables 
net present value, $ 
iteration number 


rate of utilization of auxiliary energy, Btu/hr 


rate of energy transfer at the heat exchanger, 
Btu/hr 


rate of energy transfer from the load to the 
environment, Btu/hr 


rate of useful energy collection, Btu/hr 
collector tilt correction factor 

ratio of monthly diffuse insolation on a tilted 
surface to the monthly direct insolation ona 
horizontal surface 

average ratio of monthly direct insolation on 
a tilted surface to the monthly direct insolation 
on a horizontal surface 

ratio of monthly reflected radiation on a 
tilted surface to the monthly total radiation 
on a horizontal surface 

collector tilt angle, deg 

search direction 

ambient climatic temperature, F 

collector inlet temperature, F 

collector outlet temperature, F 

storage tank outlet temperature, F 

collector loss coefficient, Btu/hr et? F 
collector loop flow velocity, ft/hr 

storage loop flow velocity, ft/hr 

lower side constraint for i-th design variable 


upper side constraint for i-th design variable 


rate of work energy used to maintain system 
flow, Btu/hr 


vector of deisgn variables 


Greek Letter Symbols 


a - absorptance of the collector absorber surface 

a* - move parameter in optimization problem 

§ - declination angle of the earth, deg 

€ - heat exchanger effectiveness 

oy - lst collector flow parameter 

bo - 2nd collector flow parameter 

n - collector efficiency 

Ge - unique flow factor (optimum F") 

c** - flow rate proportionality constant, gpm hr F/Btu 
re - ground reflectivity 


- transmittance of collector covers 


> - latitude angle, deg 
We - sunrise hour angle for horizontal surface, deg 
i We - sunrise hour angle for tilted surface, deg 


I. INTRODUCTION 


A. BACKGROUND 

The control and utilization of energy has become the 
major issue of this decade: the energy crises. The infla- 
tion rate of energy costs and particularly fossil fuel costs 
has the spector of uncontrollability as this decade comes to 
a close. Reference [1] reports a 30% increase in the price 
of oil in 1978 alone. Reports from the media in early Septem- 
ber 1979 indicate the price of home heating oil has incresed 


70% since the 1978 heating season. The crisis center appears 


to be an emerging realization that the supply of available 


energy is in fact exhaustible. This realization has spread 


from a few to the masses; unquestionably the fruit of the 


Organization of Petroleum Exporting Countries. 


In June of 1979, President Carter proposed to the Congress 


an energy strategy whereby solar energy would be providing 


20% of the Nation's energy by the year 2000. For the pur- 


poses of domestic hot water (DHW) and residential space 


heating this goal is technically achievable. The basic solar 


energy technical theory is well documented by Duffie and 


Beckman [2] and Kreith and Kreider [3] and the technology 


is continually being updated in the Solar Energy Journal [4]. 


The increased effort in solar energy research work is in 


clear evidence as the documentation in this journal appears 


to be on an exponential rise. 


The treatment of solar energy system performance is 
typically a determination of the fraction of a given heat 
load which is provided by solar energy. This fraction, f, 
is a nonlinear function of many system variables but typically 
increases in a monotonic fashion as a function of collector 
area. To simply seek maximum f and maximum collection of 
energy in an unconstrained fashion can result in extreme 
initial system costs. Accordingly, in order to assure a 
proper balance, technical performance and cost should be 
considered throughout the design process and economics should 
become an intimate part of the technical design problem; 
not an after thought at the end of system design. 

The economic ingredient in the problem of solar energy 
system design has received little attention in the literature. 
In view of this there is little information in the literature 
on how to proceed with design improvement of a solar energy 
system. Accordingly, a broad objective for this work was 
to develop a solar energy system model including economic 
considerations and to seek system design improvement by 
using the resulting eco-technical model together with an 


optimization algorithm. 


B. SCOPE 

The conversion of solar energy into useful work covers 
a broad spectrum of collection schemes. These include thermal, 
photovoltaic, biological, wind, and ocean thermal energy 


conversion. Reference [5] is an excellent introductory 


information source covering many methods of solar energy 
conversion and has become a classic. This thesis effort 
was limited to the thermal conversion of solar radiation 

for the purpose of generating low temperature heat. Since 
low temperatures are involved, only flat plate collectors 
were studied. Flat plat collectors can provide temperatures 
on the order of 200° F compared to concentrating collectors 
which can provide temperatures as high as 6700° F. The 


restrictions on temperature were chosen in keeping with 


availability analysis; matching the solar conversion system 


to its task. The task of providing high entropy, low quality 
DHW and space heating energy is matched with low availability 
(i.e., minimum utilization of available energy in high 
quality, low entropy forms such as oil, gas, and central 
power). The collector working fluid was restricted to water. 
Further, since the collection of solar energy is a stochas- 
tic process, Only analysis for long term performance was 
studied. The analysis of the dynamic performance of specific 
systems in response to hourly climatic data has previously 
been conducted by Kline et al [6] and now forms the basis 


for long term analysis. 


C. OBJECTIVE 

The objective of this thesis was to develop methodolo- 
gies for: 

l. System synthesis including a determination of optimum 


system design variables including collector area, collector 


tilt angle, heat exchanger sizes, and flow rates. This 
analysis would be based on using shelf item collectors 
characterized by performance parameters FL.to and FUL 
which have been determined by tests conducted in accordance 
with refs. [7 and 8]. 

2. Collector design optimization including a determin- 
ation of the optimum combination of collector design 
variables such as materials and geometries; and which 
would yield the optimum system performance. 

Fundamental to the effort was the utilization of: 

1. COPES/CONMIN (Control Program for Engineering 
Synthesis with Constrained Function Minimization), a design 
improvement algorithm developed by Vanderplaats [9, 10, ll 
and 12]. 

2. Standardized long term solar energy load fraction 
correslations (f - CHART) developed by Kline [6]. 

3. Economic considerations combined with technical 
analysis. 


4. Long term climatic data obtained from (NOAA) 


National Oceanic and Atmospheric Administration [13]. 
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II. SYSTEM MODEL 


A. OVERALL SYSTEM DESCRIPTION 

A schematic diagram of the solar heating system is 
shown in Figure l. The system consists of four heat 
exchangers, a storage tank and associated system piping. 
The first heat exchanger is called the collector and is 
used to transform incident solar radiation into useful 
thermal energy. At the second heat exchanger the collected 
energy is transferred to the secondary loop and stored in 
the form of sensible energy in the storage tank. The 
sensible energy is transferred from storage to the loads 
via the third and fourth heat exchangers. Energy transfer 
is achieved via four forced convection flow loops. fMThis 
model was further simplified by assuming total utilization 


of the energy stored in the tank and thereby eliminating 


the need to detail the DHW and space heating heat exchangers. 


B. ENERGY BALANCE 


The steady state power balance for the system is 


Q. + WwW (1) 


5 + aye a 


L 


where, 


Q. 2 Useful rate of energy collection (Btu/hr) 


Ww = Work energy utilized to maintain system flow 
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Q = Rate of utilization of auxiliary energy 
which is required in addition to solar 
energy in order to meet the DHW and space 
heating load 


Qy = Rate of energy loss by the system 


For the purposes of this thesis, it has been assumed 
that Wp is a known or estimated parameter. Most analyses 
in the literature do not account for W_ since it is typically 
assumed to be a very small fraction of the energy required. 
The percentage of the energy load which is supplied 


by solar energy is defined as: 


£ = 2 (2) 


where the summation is taken over 12 monthly values. 

£ is a monotonically increasing function of collector 
area. No optimum f is evident. Figure 2 depicts a typical 
f distribution resulting from a design problem. Since 
initial costs are strongly dependent on collector area (see 
Section III), attempts to seek maximum f can result in very 
high initial system costs. At this point the motivation to 
include economic considerations as a part of the system 


design model become clear. 


C. ECONOMIC ANALYSIS 
The matter of economics has been treated as it evolved 
in this study: a simple common sense approach. This 


approach includes: 


l. An economic equation which utilizes technical 
parameters obtained from mass, momentum and energy balance 
considerations together with economic parameters. 

2. Standardization of the economic equation by char- 
acterization of all terms in present dollar values. 


The economic analysis is similar to an energy balance. 


* 
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Considering the economic system as the monetary account 
of a consumer and the time period as the economic lifetime 


of the energy conversion equipment: 


Pout - D storage 


The value if dollars of the energy 
produced by the equipment 


The expense dollars to obtain, operate 
and maintain the equipment 


Dstorage = The net savings 


The rational consumer seeks maximum D This 


storage’ 
methodology is a simple capital budgeting technique and is 


contained in any standard reference in accounting or finance. 


D. THE OBJECTIVE FUNCTION 

Implicit in the use of the economic model is an assump- 
tion that the terms in the equation can be expressed as 
variables of the technical functions which have been 


defined by the physics of the engineering problem. For 


Bh ew en ebihts sarorw ind oie bloom = ~~ 


the solar energy system utilized for domestic hot water 


and space heating the transformation is simple: 


Net Savings = Fuel Savings - System Costs 


This equation becomes the objective function by 


standardizing all dollar amounts to present values. The 


a a i nl a ee 


resulting equation is: 


NPV = fo,C.F} ~ A.C, (4) 


where, 


NPV = Net present value of the solar investment, $ 


Cc = Fuel cost for the energy replaced by solar 
energy, $/Btu 


Fi = Present worth factor which standardizes fuel 
savings during system life into present 
dollars, years 


. 2 
A, = Collector area, ft. 
Cy = System Lifgcycle cost per unit collector 
area, $/fts 


The effort now proceeds to determining each of the 
above dependent parameters and to obtaining those combina- 
tions of these variables which will result in maximum 


values of NPV. 
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III. SYSTEM OF EQUATIONS 


A. GENERAL 

The objective function has been fully developed by 
examining each of its parameters individually. Completion 
of this development has resulted in a non linear function 
of eight independent design variables; a hypersurface in 
eight space. When formulated in this manner, the problem 
is clearly not amenable to analytical solution. Ina very 
simple manner however the problem can be interfaced with 
the design improvement algorithm COPES/CONMIN (see Section IV). 
The primary contribution for this effort comes from Kline [6] 
who developed the f-chart correlations for determining monthly 
sOlar energy load fractions fj: The annual load fraction 
fF is then computed by a weighted average of the monthly 


values. 


B. SPECIFIC PARAMETERS 
1. Annual Energy. Load, Qy 

The annual energy load consists of two basic ingredi- 
ents; the DHW load and the space heating load. Assuming an 
average mean ground temperature the month to month DHW load 
was considered constant except for the variation due to month 
length. The methods for determining the space heating loads 
typically follow the ASHRAE manual, reference (l4]. Refer- 


ence [(15] is recommended as a text for space heating load 


ts 


computations. For the purposes of this study representative 
building loss parameters were used. Typical building con- 
ductance or (UA). values of 30,000, 20,000 and 10,000 
(Btu/deg F day) have been used based on a standard building 
of 1750 ft? of floor area, a heat transfer surface of 5000 
ft? and building loss coefficients, UL. Of .25, «17; and «09 
Btu/ft* hr F respectively. The conductance, UA, is the 
space heating load at design conditions which has been 
estimated in the manner of Ref. [14] and divided by the 
design temperature difference. 
2. Annual Solar Energy Load Fraction, = 

fF is a function of collector performance parameters, 
thermal physical properties of the working fluids, fluid 
flow rates, heat exchanger performance parameters, collector 
tilt angle, climatic conditions and latitude. The climatic 
conditions determine the load distribution (heating degree 
days), the solar energy flux distribution (insolation) and 
the environmental stress on the collector (ambient tempera- 
ture). 

a. Collector Performance 

A determination of collector performance proceeds 

from the well known collector equation of Hottel and Whillier 
{16) which is the result of an energy balance on the collector: 


we = FLUPAS (IRta = Up (tT, = T,)) (5) 


where, 
Useful rate of energy collection (Btu/hr) 


Collector heat removal factor (dimensionless) 


2 


Collector loss coefficient (Btu/hr ft" F) 


Collector ares (ee?) 


. 
Solar insolation on a horizontal surface (Btu/hr ft”) 
Collector tilt correction factor (dimensionless) 
Transmittance of collector cover system (dimensionless) 


Absorptivity of the collector absorber plate 
(dimensionless) 


Collector fluid inlet temperature (F) 
Climatic ambient temperature (F) 


The analytical expressions for e. and U, are complex 


L 
and result from a lengthy development which can be found in 


references [2 and 3]. The computation of v. and U, using 


L 
these expressions is not necessary for system design analysis 


Se ee 


using shelf item collectors since the parameters FOL and 


Pita can be obtained from collector efficiency tests. 


Collector efficiency is defined as: 


Pe TN RE ae eB Ee 


Assuming that U, is not sensitive to changes 


L 
in the environment and combining equations (4) and (5) results 


n 2 Puta - FU, ((T, - T,? / IR) 


Equation (7) is linear in the collector 


parameter (T; - T,)7IR and forms the basis for determining 
collector performance. It is clear that the parameter 
FU, is obtained from the slope and Fit is obtained from 
the intercept of equation (6). Some typical curves of 
collector efficiency are included in Figure 3. 

Two standards have been developed for collector 
testing to measure collector performance, namely references 
[7] and [8]. The standards differ in the independent 
variables used in the performance characterization. 
Reference [7] the ASHRAE standard uses T = T, = Ti 
ref (8] the NBS standard uses T = (T; + Ty) 72, where T, 
is the collector outlet temperature. All collector 
performance test data used in this effort is based on the 
ASHRAE standard. All collector data used in this effort 
is based on actual collector tests conducted by NAVFAC 
and reported via ref [17]. 

b. Collector to Storage Coupling 

The collector becomes coupled to the storage 
tank via the collector to storage heat exch ‘er following 
the method of de Winter [20]. It is assumed that the rate 
of energy transfer at the heat exchanger is equal to the 


rate of useful energy collection: 


Qhex = Q. = e¢(mc_) 


where, 
= = Maximum system temperature or collector outlet 


7 = Minimum system temperature or storage outlet. 


Equations (5) and (8) are combined to eliminate 


T, and T, dependence to give: 


Qy = PLU, (IRTa - UP (TS - T,)) (9) 
where, 
rS Pex ¢ (10) 
f FLULAQ Ce \ fe 
(me,) . € (MC) min ; 


c. Solar Insolation at Optimum Tilt Angle, Ty 
The solar radiation intensity on a tilted 


surface is by definition: 


&, * IR (12) 


Solar insolation is typically measured and 
reported for a horizontal surface. Improvements in a 
solar collector installation are sought by tilting the 


collector to the optimum angle. The computation for 
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optimum angle is iterative for each design situation. 
Several rules of thumb have evolved: 
(1) Collectors should be oriented at a slope 
of .9 times the latitude angle for maximum annual collection. 
(2) For DHW systems where loads vary little 
during the year the best angle of tilt is equal to the 
latitude angle. 
(3) The optimum collector tilt angle for least 
; cost per Btu delivered for building heating is approximately 
the latitude angle plus 15 degrees. 

The methodology for determining collector tilt 
angle used in this effort was as follows: 

(1) Horizontal monthly data was obtained for 
97 locations from NOAA in reference [13]. 

(2) An algorithm was developed to compute 
monthly horizontal extraterrestrial radiation intensities. 
This computation is a function of latitude, daily hour 
angles on horizontal and tilted surfaces, daily declination 
angle, and collector tilt angle. The azimuth angle was 
always chosen for a due south collector orientation. 

(3) The methodology follows that of Liu and 
Jordan in reference [19]: 


R = (l1- QR + R, + R (13) 


9 
¢ = 1.3903 - 4.0273 (4) + 5.5315(4)7 - 3.108 (+) ° (14) 
(eo) ° (eo) 
w_'sin(g-s)sin § + cos($-s)cos § sin w_' 
i ww ut : (15) 
da w Sin > sin § + cos § cos § sin Ws 


w ' = min{cos’-(-tan(¢-s) tan o}y cos + (-tan d tan 6)] (16) 


Ws * cos }[{-tan » tan 8] (17) 
i * = (1 + cos s) (18) 
1 
R = #% (1 = cos s)o (19) 
where, 
R = monthly slope correction factor 


Ry = average ratio of monthly direct radiation on a 
tilted surface to monthly direct radiation on a 
horizontal surface 


Ra = ratio of monthly diffuse radiation on a tilted 
surface to monthly diffuse radiation on a 
horizontal surface 


R. = ratio of monthly reflected radiation on a tilted 
surface to monthly total radiation on a horizontal 


surface 

3 = ratio of monthly diffuse radiation on a horizontal 
surface to monthly total radiation on a horizontal 
surface 

-* ratio of monthly total radiation on a horizontal 

fo) surface to monthly total extraterrestrial radiation 


On a horizontal surface 
Ww, ™ sunrise hour angle on a horizontal surface 


w,'= sunrise hour angle on a tilted surface 
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s = collector tilt angle 


> = latitude angle at collector site 
§ = declination angle of the earth 


© = ground reflectance for area adjacent to collector. 


3. Fuel Cost Parameters 
Fuel cost parameters included a selection from among 
3 fuel types (i.e. oil, electricity, or gas), unit of issue 
cost, fuel heating value and efficiency of the auxiliary 
heating system. The SOLOAD-1 system allows complete freedom 
in the selection of these parameters. A typical set of fuel 


cost parameters used in a design experiment included: 


Fuel: Oil 

Unit Cost: .9 ($/Gal) 

Heating value: 142,000 (Btu/Gal) 
Efficiency: "i | 

Resulting Fuel Cost: Ce = $9.05 per 10° Btu 


4. System Cost Parameters 


System cost considerations typically include 
initial costs and life cycle costs for operation and 
maintenance. Detailed guidance for complete cost consi- 
derations is included in ref. [3]. Further, reference [20] 
contains guidance to be used in determining costs for US 
Navy installations. The following cost parameters were 


used in SOLOAD-1: 


Total system life cycle cost 


Total cost per square foot 


Me Ae OEY RM EEE othe 


Initial system cost 


r- 


Operation and maintenance costs 


< 


Collector cost 


a 


System installation cost 


Ht 


Storage tank costs 


cr 
ro 
| 


Collector to storage heat exchanger costs 


“a 


Cc 
c 
Cc 
Cc 
Cc 
C 
Cc 


Pumping power costs 


ue) 


Maintenance cost factor 


a | 
3 
vy 


" 
ro 


Present worth factor 


>, ,7N 
. bs tk+a9 
rs 


where, 


So CRE 


et ep BS (25) 


Pee MT Me 3 


i = Discount rate 


oa Nara 


Fuel inflation rate 


N = System life in years 


It should be noted that Fi occurs explicitly in 
the first term of the objective function and implicitly 
in the second term. It should also be noted that initial 


costs are not amortized in the SOLOAD-1 model. 
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IV. NUMERICAL OPTIMIZATION 


A. SIMPLE OPTIMIZATION 

Design problems typically seek the minimization or 
maximization of an appropriate parameter within a frame- 
work of constraint specifications. The parameter to be 
optimized may be a function of several design variables and 
is termed the objective function. Other parameters which 
may be separate functions of design variables must not 
exceed specified bounds for the design to be acceptable. 
These parameters are termed design constraints and are not 
to be confused with limits which may be set on design 
variables which are usually termed side constraints. 

Engineering problems can be numerically coded for an 
analysis (once through) solution. The simplest scheme for 
Optimization may consist of a series of loops through the 
computer code which may cycle many combinations of design 
variables. The combination of variables which provides the 
best design and which also satisfies the constraints is 
then considered optimum. This approach may be economical 
for a design problem with just a few design variables and 
short computer time requirements. A design problem with 3 
design variables, ten values for each design variable, and 
one-tenth seconds central processing unit (CPU) time for 


each analysis would take a total of 100 seconds of CPU time. 


The solar energy optimization design problem as characterized 


by the code developed by this thesis has a minimum of 8 
design variables; each analysis of its objective function 
requires about 2 CPU seconds. Using the simple approach 
and assuming ten values for each design variable would 
result in a CPU time of 68 years for each design problem. 
Clearly, a more rational approach to optimization is 
necessary. 

Vanderplaats [9] suggested that many special algorithms 
for numerical optimization have been proposed in recent 
years, but that in many cases unsuspecting practitioners 
find their particular optimization problem unsolved only 
after large amounts of time and effort are expended. This 

: can occur usually because of inexperience by the practitioner 
in determining the limitations of specified algorithms. 
Vanderplaats (10, ll, and 12] has developed a FORTRAN 

coded algorithm capable of optimizing a very wide class of 
engineering problems. The system includes COPES (Control 
Program for Engineering Synthesis) and CONMIN (Constrained 
function Minimization). This optimization system is referred 


to as COPES/CONMIN. 


B. THE COPES/CONMIN SYSTEM 
1. Terminology 
CONMIN is a FORTRAN program in subroutine form for 
the solution of linear or non-linear constrainted optimiza- 


tion problems. The user prepares an analysis program. The 
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program must be named SUBROUTINE ANALIZ. The process of 


computer aided design or of trade off studies with a minimum 
of man-machine interaction becomes fully automated via the 
COPES program. Three basic definitions are required: 

Design Variables. Those parameters which the 
Optimization program (CONMIN) is allowed to change in order 
to improve the design. Design variables appear only on the 
right hand side of equations in the analysis program (ANALIZ). 
Limits imposed on design variables are termed side constraints. 

Objective Function. Usually the single parameter 
which is to be minimized or maximized during optimization. 
The objective function always occurs on the left side of 
the equation in the analysis program. (Refer to [12] for 
exceptions.) The equation defining the objective function 
may be linear or non-linear, implicit or explicit, but must 
be a function of the design variables to be meaningful. 

Design Constraints. Any parameter which must not 
exceed specified bounds for the design to be acceptable. 
Constraint parameters always appear on the left side of the 
constraint function equations. Constraint functions may be 
linear or non-linear, implicit or explicit, but must be 
functions of the design variables. 

Assuming that the optimization process requires the 
maximization of a particular objective function, the general 


Optimization problem can be stated as: 
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a. Find: X which maximizes f (X) 
b. Subject to: 
(1) CONSTRAINT EQUATIONS G; (X) < 0, j = 1, NCON 
(2) SIDE CONSTRAINTS VLB < Xs < VUB., i=1, NDV 
Where, X = X(X), Xo, +++, X1) is the vector 


of NDV design variables, F(X) is the objective function and 


G(X) are the set of NCON constraints. VLB, and VUB, are 


variable. 


2. Methodology 


The solution process proceeds as follows: 


: 
' 
§ 
the lower and upper bounds respectively on each design 
| 
: 


a. The user prepares an analysis subroutine which ; 
defines X, F(X) and G(X). This subroutine must be named 
ANALIZ. ANALIZ must have three segments; input, analysis 
ana output keyed to COPES flags, ICALC = l, 2, or 3 
respectively. 

b. The user prepares an input data file for COPES 
which includes a wide variety of system options, appro- 
priate matching mechanisms between ANALIZ and CONMIN and 
the constraint boundaries. 

c. Using the initial vector of design variables 
COPES obtains an initial solution from ANALIZ and subsequent 
solutions by updating X as determined by CONMIN. Any analy- 

- sis sOlution which satisfies the constraint equations and 


the side constraints is a feasible design. If an analyti- 
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cal solution violates any of these constraints it is an 


infeasible design. The minimum feasible design is optimal. 
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The feasibility determination includes: 

(1) If a constraint equation is violated 
tives, 12 G(X) > 0) then the jth constraint is violated. 

(2) If a constraint equation equality 
condition is met (i.e., if G, (X) = 0) then the jth con- 
straint is active. 

(3) If a constraint equation condition 
is met (i.e., if G5 (X) < 0) then the constraint is inactive. 

Note that CONMIN is designed to minimize objec- 
tive functions; the process of maximizing an objective 
function is concerned with minimizing the negative of an 
objective function. 

(4) Similarly, side constraints may be 
inactive or active but side constraints will nevcr be 
violated in a particular analysis computation because 
they are specified limits not dependent variables which 
is the case for design constraints. 

(5) All inequality conditions are repre- 
sented by a band around the zero condition due to computer 
limitations in defining zero. 

d. If the initial analysis solution is infeasible 
CONMIN moves towards a feasible solution by adjusting the 
design variables appropriately. The optimization process 
then proceeds in an iterative fashion as follows: 


(1) The iterations are governed by the 


recursion relation: 


eg 


ee Aree cm ree 


where, 


q = iteration number 

a* = a scalar move parameter which defines the 
distance of travel in the direction of 
search 


direction of search 


wl 
" 


(2) § is determined such that F(X) will 
be minimized without violating any constraints. CONMIN 
calculates the gradient of F(X) by using finite difference 
techniques. Because no constraints are violated, the 
greatest improvement in F(X) will be realized by moving 
in the direction of steepest ascent, the gradient of F(X). 

(3) Once the directive is known, the 
move parameter, a*, which will allow the largest magnitude 
improvement in F(X) is to be found. A one dimensioned 
search of the § direction is carried out until the best 
improvement number, a*, is found. 

(4) CONMIN utilizes methods other than 
the method of steepest descent in determining S$ particularly 
in the presence of active constraints. These methods 
include the method of conjugate directions developed by 
Fletcher and Reeves [21] and the method of feasible direc- 
tions developed by Zoutendijh [22] and implemented by 


Vanderplaats and Moses [23]. 
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e. CONMIN continues to iterate for an optimal 
design by computing successive § and a* values always 
keeping within the defined constraints. If there is no 
relative or absolute change in F(X) for three successive 


iterations, the optimum is considered found. 


Soa aaeieednnenaaietadaieadi di seeesaendaaitaamneeememenina ae 
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V.  SOLOAD-1 SYSTEM 


A FORTRAN coded algorithm for the analysis of systems 
containing shelf item collectors has been developed. The 
algorithm has been named SOLAR ENERGY OPTIMIZATION ANALYSIS 
OR DESIGN (SOLOAD-1). The portions of the algorithm 
reserved for design optimization of collectors were not 
completed and have not been included in SOLOAD-l. The 
interface between SOLOAD-1 and COPES/CONMIN has been 
developed via the Control Program (CP-67) and the Cambridge 
Monitor System (CMS) a time-sharing system developed for 
the IBM System 360 Model 67. SOLOAD-1 includes: 

A. Three executive programs to provide the interfaces 
among the system elements and to initiate design problems. 

B. An analysis subroutine which includes the system 
of equations covered in Section III. A summary of the 
objective function, design variables, and design constraints 
is included in Appendix A. 

C. Three auxiliary subroutines including a calendar 
data array, climatic data and user defined input parameters. 
D. Data files for user input to SOLOAD-1 and COPES/ 

CONMIN. 

The interface among SOLOAD-1 system elements and 
COPES/CONMIN is shown in Figure 3. SOLOAD-1l programs and 
subroutines are included in Appendix B. Data files are 


included in Appendix C. 


ee 
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VI. RESULTS 


A. GENERAL 

The first objective of developing a methodology for 
the analysis of systems with previously defined collector 
design was achieved. The second objective of developing 
a methodology for optimizing the collector component design 
was not achieved. However, it appears that the potential 
exists for using the system analysis methodology to obtain 
improvements in collector design. For example, if cost 
data on various types of collector cover plates or collector 
absorber surfaces can be correlated with collector per- 
formance parameters, then the system model could be used 


for quick checks on collector design improvements. 


B. VERIFICATION 

A single design problem was selected and started using 
10 different sets of the starting design variable vector. 
The same optimal design result within + 0.5% was achieved 
for each run. This was the only actual verification deemed 
possible at this time due to lack of any known optimization 


data in the literature. 


C. EXPERIMENTS 
Approximately 50 design experiments were analyzed. 
Each problem was characterized by a unique identification 


number which could tie together the location, economic 
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environment, the collector, and the space heat transfer 
coefficient, UA, for the experiment. For each experiment 
an input parameter summary report and output parameter and 
results summary report was generated. These reports are 


included in Appendix D. 


D. SPECIFIC RESULTS 
1. Constraints 
The only active design constraint in most experi- 
ments was the tube thickness in the collector loop at the 
heat exchanger. Some experiments resulted in no active 
design constraints. There were no active side constraints 
in any experiment. 
2. Collector Flow Factor 
The collector flow factor F" is defined as the 
ratio of the collector heat removal factor to the collector 


efficiency factor: 


The flow factor result obtained in all optimization experi- 
ments was the same. F" = .948 + .008. 
3. Collector Tilt Angle 
Each experiment conducted for a particular geo- 
graphical location resulted in the same collector tilt 
angle. However, experiments for different locations but 
the same latitude resulted in significantly different 


optimum collector tilt angles. 


39 


VII. CONCLUSIONS 


A. The SOLOAD-1 system in conjunction with COPES/CONMIN 


appears to offer the potential for further improvement and 


Ee 
ee 
—— 


potentially a valuable automated technique for solar energy 
system design. 

B. The "rule of thumb" typically used for collector tilt 
angle optimization (i.e., latitude plus 15°) should be 
used with caution since preliminary results indicate a 
strong tendency for climatic dependency. 

C. The uniqueness of the flow factor was suggested by 
the continued result of .948 + .008 for each experiment. 
Pending confirmation by further testing it appears that 

a simple correlation for determining optimum collector 
loop flow rates may be available. Based on a flow factor 
of 0.948 the resulting correlation would be: 


G = .01955 r— (gpm/£t? of collector area) 


L 


This correlation is developed in Appendix E. 
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VIII. RECOMMENDATIONS 


A. The basic model of SOLOAD-1 should be upgraded to 
include specific characterization instead of simple 
parameter selection and input for: 

1. Pumping power in all four loops 

2. Inclusion of the DHW heat exchanger 

3. Inclusion of the space heat exchanger 

4. Inclusion of building loss coefficients as a 

design variable. 

B. Prior to any additional experiments with the present 
model, a complete survey of the industry should be con- 
ducted for collector parameters and costs. Following this 
survey, a series of experiments should be conducted in 
search of correaltions among collector performance parameters 
(heGe, PL (ta) and FUL), collector unit cOsts, and system 
performance. 
C. A larger sample size should be used to verify the 
uniqueness of the optimum collector flow factor (F") as 
suggested by the results of 50 experiments. 
D. The remainder of the NOAA climatic data bank (i.e., 67 
more cities) should be included in SOLOAD-l. The optimum 
collector tilt angle could then be computed for the 97 
NOAA regions. 
E. The model should be upgraded to accommodate analysis 


for systems with air as the working fluid. 


F. The model should be upgraded to accommodate amorti- 
zation of the installation cost instead of just initial 


cash payment. 
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FIGURE 3 Typical Collector Efficiency Curves 
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APPENDIX A 


SUBROUTINE ANALIZ SUMMARY 
OBJECTIVE FUNCTION 
— + a 
NPV = f Q, Ce Fi A, Cs 


DESIGN VARIABLES 


Collector Area Aa 
Collector Tilt Angle s 

Collector Loop I.D. a; 
Collector Loop O.D. a, 
Heat Exchanger I.D. aoa 


Collector Flow Velocity v 
Storage Flow Velocity Vv 


Heat Exchanger Length L 


DESIGN CONSTRAINTS 


= dy5 = dy 


= Reynolds number, Collector loop 


= Capacity ratio, Coren 


1 
2 
3 

G, = Reynolds number, Storage loop 
3 
, * lst Flow parameter, i he, /AF 
- 


= 2nd Flow parameter, bo * the ,/AF 


t 
GU, 
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APPENDIX B 


SOLOAD-1 Computer Program 
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APPENDIX D 


EXPERIMENT REPORT SUMMARIES 


Each optimization problem or experiment is represented 


an a i sel asi ll 


by an input summary report and an output summary report. 
These reports are tied together by a unique identification 


as follows: 


6 SIE RE EDN ie rina: 


elle pe ee 
| 
where, ¢ 
: 
LL = location identification number in 4 
accordance with Appendix C : 
N, = present worth factor identification number ; 
1- Fi = 18.22 ; 
2-P! = 24.34 i 
i i 
ft 


N = collector identification number 
1 - Solarnetics 


2 - American Sun 


es 


3 - Federal Prison Institute Double Glaze 


N, = heat load conductance 


1 - VA = 30000 Btu/HDD 


CTE er Beg RO HNO IE 


2 - UA = 20000 Btu/HDD 


3 - UA = 10000 Btu/HDD 


FEO Ee a aS YN Ee. 


76 


The following report sets are included: 


1111 
L112 
1113 
1213 
1223 
1232 
1233 


SLi 
9112 
9113 
9213 
9221 
9222 
9223 
9231 
9232 


LSlil 
13112 
13223 
13232 


2222 
2223 
2232 
2233 


20111 
10112 
10112 
10213 
10221 
10222 
10223 
10231 
10232 
10233 


14111 
14112 
14224 
14232 


S31ii 
3112 
3213 
3222 
3223 
3232 
32323 


ALiL1L 
Lii12 
LLLis 
11221 
L1Z22 
11223 
LL232 


LoLii 
LS112 
iS2Z3 
15232 


4111 
4112 
4113 
4222 
4223 
4232 
4233 


LZiil 
12112 
L4ii3 
L2221 
12222 
12223 
12232 


16111 
16112 
16222 
16223 
16232 
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APPENDIX E 


POTENTIAL CORRELATION FOR OPTIMUM COLLECTOR FLOW RATE 


Reference [2] defines the collector heat removal factor, 


F as: 
zr * 


-F'U 
Ge 
P 


Ge 
FL = GP i - exp( dy} (Bl) 


Further, the collector flow factor, F", is defined as, 


¥. Ge “Fu, 
) Fr = are {l - exp( Se (B2) 
L P 


Reference [2] shows that F" approaches the unity value 
asymptotically as the parameter Ge, /U,F* increases. 
Reference [3] recommends capacity rates of 10-15 lb/hr fe 
(.0223 = .0334 Gpm/£t?) as the best compromise among 
collector heat transfer coefficient, fluid pressure drop, 
and energy delivery. The unit ft. refers to collector area. 

The results of this thesis suggested a correlation 


between the collector performance parameter FU and the 


L 
capacity rate. Once these results are verified by further 
testing including model changes to include fluid pressure 
drop, a simple correlation may be available and follows 


from: 
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Let, x* = The optimum flow factor determined from 
computer experiments. 


r, = Ge_/F'U 
c,/ 


*1 L 
t= Go, /F Uy 


It follows from equations (Bl) and (B2) that: 


me 
2 


— =i4t.l- exp(-1/;,) 
or, 


hate l 
» = Fi = 
1 In(1l 1/35) 


x 
Let, $$) = The optimum parameter which corresponds to 
<* and 1S obtained by solving equation (B2) 


for 1 
Therefore, 
fh 
<* = lL 
52 
or 
sf 
$2 - Pz 
or 


8 
Me] 

Ne 

i’ | 
n" 

Gc 
0 

it 
"1 


Based on the results of a limited number of computer 


experiments all yielding an apparent invarient flow factor, 
F" = 1948: 


x** = .01955 (Gpm hr F/Btu) 


Or, 


2 
G = -O1955F_U, (Gpm/ft-) 
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